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ABSTRACT 

Wireless OFDM channels can be approximated by a time varying fil- 
ter with sparse time domain taps. Recent achievements in sparse sig- 
nal processing such as compressed sensing have facilitated the use 
of sparsity in estimation, which improves the performance signifi- 
cantly. The problem of these sparse-based methods is the need for 
a stable transformation matrix which is not fulfilled in the current 
transmission setups. To assist the analog filtering at the receiver, 
the transmitter leaves some of the subcarriers at both edges of the 
bandwidth unused which results in an ill-conditioned DFT subma- 
trix. To overcome this difficulty we propose Adaptive Threshold- 
ing for Sparse Signal Detection (ATSSD). Simulation results confirm 
that the proposed method works well in time-invariant and specially 
time-varying channels where other methods may not work as well. 

Index Terms — Sparse channel. Adaptive thresholding, OFDM 
channel estimation, MMSE 

1. INTRODUCTION 

Recent standards such as [?, ?, ?] introduce OFDM communica- 
tion as one of the best options for wireless transmission of multime- 
dia signals including video. OFDM transmission, although resistant 
against multipath fading, requires accurate estimation of the Channel 
Frquency Response (CFR) at the receiver for appropriate decoding 
of the data. 

For mobile reception, the transmitter reserves some of the sub- 
carriers in each OFDM symbol (which may vary in different sym- 
bols) for a predefined pattern of data, referred to as Pilot tones, 
which are a priory known at the receiver. Since in OFDM transmis- 
sion the data are directly sent in the frequency domain, these pilot 
tones provide the receiver with irregular noisy samples of the CFR 
in the time interval of that symbol. The main task of the channel 
estimator block of the receiver is to estimate CFR at the non-pilot 
subcarriers using the obtained samples at pilot tones. Conventional 
methods involve some sort of interpolation between the obtained 
samples, ranging from the simple Linear Interpolation (LI) to more 
complicated approaches like splines. Under limited time spread and 
Doppler frequency, the time varying OFDM channel can be viewed 
as a lowpass 2D signal [?]. In this sense the problem of channel es- 
timation becomes equivalent to the reconstruction of a 2-D lowpass 
signal from its irregular samples [?]. 

Recent results in sparse signal processing such as Compressed 
Sensing [?], paved the way to exploit the inherent time domain spar- 
sity of the channel in its estimation [?, ?, ?]. Using l\ minimization 
algorithms (Basis Pursuit) introduced in [?], reduction of the num- 
ber of required pilots for channel estimation is proposed in [?] while 



in [?] a sparsity-based estimation method employing the Matching 
Pursuit is devised. 

A main drawback of these methods is that they do not consider 
zero padding in their scenario. In current OFDM standards the band- 
width is not fully occupied; a number of the subcarriers at both edges 
of the bandwidth are set to zero (hence the name zero padding) to in- 
crease the allowable transition band of the analog bandpass filter at 
the receiver Zero padding results in an unstable frequency to time 
transformation (ill-conditioned transformation matrix). In addition, 
due to the lack of pilot in zero padded parts, common time domain 
techniques are impractical. In this paper we propose a sparse chan- 
nel estimation method called Adaptive Thresholding for Detection of 
Sparse Signals (ATSSD) that works adequately even for zero-padded 
OFDM systems. Simulation results using the DVB-H standard con- 
firm that ATSSD estimates almost the exact CFR in time-invariant 
channels. Furthermore, unlike the previous methods, the perfor- 
mance under varying channels degrades only slightly as the Doppler 
frequency increases. 

2. SPARSE OFDM CHANNEL ESTIMATION 

To show how the aforementioned challenges in OFDM channel es- 
timation can be overcome using the inherent channel properties, we 
will briefly explain the problem in two subsections. In 12. II OFDM 
system components related to channel estimation are reviewed. In 
12.21 in addition to restating the OFDM channel estimation as a sparse 
signal processing problem, we will explain the difficulty of partial 
use of bandwidth (zero padding) in estimating sparse channels. 

2.1. Review of OFDM System Components 

At the OFDM transmitter, coded binary data are interleaved, grouped 
and mapped to constellation points based on a specific modulation 
scheme. The resulting complex symbols which will form the fre- 
quency domain of the transmitting signal are rearranged into blocks 
and for each block a number of pilots are inserted among the data at 
predefined locations. To facilitate the RF analog filtering at the re- 
ceiver, these blocks are zero padded to form the final blocks (OFDM 
symbols, denoted by X'"' where n is the block number); i.e., the 
bandwidth is not fully occupied. The resulting symbols are con- 
verted into the time domain (x^"') by the use of the IFFT algorithm, 
where each symbol is extended by a Cyclic Prefix (CP, a copy of 
the last part of each OFDM symbol) and these final blocks are se- 
rially transmitted. The cyclic prefix is inserted to prevent possible 
Inter-Symbol Interference (ISI) in existence of multipath channels. 

At the receiver, the cyclic prefix is removed and the received sig- 
nals, y(n), are sent to an FFT block. When the duration of the chan- 



nel impulse response is shorter than the cycUc prefix length, no ISI 
is expected between OFDM symbols. Assuming also no synchro- 
nization error, the equation relating the n"* transmitted and received 
symbols, respectively X(n) and Y(n), is: 

Y(n) =X(n)0H(n)+W(n) (1) 

where represents element by element multiplication of two vec- 
tors, H(n) is a vector containing the samples of the Channel Fre- 
quency Response (CFR) corresponding to the n*'' OFDM symbol 
and W denotes the sampled vector of the AWGN noise in the fre- 
quency domain. Since the receiver is priorly aware of pilot positions 
and values, it can obtain a noisy estimate of the channel frequency 
spectrum at these subcarriers. In this sense, the OFDM channel is 
similar to a 2-D lattice in the discrete time-frequency plane from 
which certain points are known due to pilots; the goal of channel 
estimation is to estimate the rest via interpolation. The interpolated 
channel estimate is then used in the equalization block to obtain an 
estimate of the transmitted constellation points (X(n)). After equal- 
ization, the approximated OFDM symbol is demodulated, deinter- 
leaved, and decoded to produce the binary output data. 

2.2. Restating OFDM Channel Estimation as a Sparse Problem 

Considering the sparse distribution of the scattering objects, the 
OFDM channel becomes sparse in the time domain. Thus, by ex- 
ploiting this sparsity, a better estimate could be obtained using the 
time domain. The resulting estimate is then transformed into the 
frequency domain by use of the FFT algorithm. In this sense, the 
problem of channel estimation becomes equivalent to finding the 
sparse Channel Impulse Response (CIR) (h) from the equation: 

Hp = Fp • h + Wp (2) 

where Hp is the vector of observed channel coefficients at pilot sub- 
carriers, Fp is the sub-matrix of the DFT matrix obtained by keep- 
ing the rows of the FFT matrix that correspond to pilot positions 
and Wp is the frequency-domain noise vector at pilot positions. As 
stated earlier, in the case of no ISI, the length of the channel cannot 
exceed the length of the cyclic prefix (Ncp). Thus, ^ can be fur- 
ther simplified in the zero ISI case; only the first Ncp elements of 
h (hcp) can have non-zero values. 

Hp = Fp,cp ■ hcp + Wp (3) 

where Fp,cp is the sub-matrix of the DFT matrix obtained by keep- 
ing only the first Ncp columns of Fp. Recently the idea of using 
time domain sparsity in OFDM channel estimation was proposed in 
[?] for decreasing the number of pilots. The authors proposed com- 
pressive sensing algorithms to find the sparsest time domain chan- 
nel. They proved that in case of uniform pilot insertion, the matrix 
Fp,cp in H) satisfies the uniform uncertainty theorem described in 
[?]. As a result, linear-programming-based algorithms used in com- 
pressed sensing, similar to the ones introduced in [?] can be applied 
to OFDM channel estimation. However, the authors of [?] did not 
consider zero-padding at the endpoints of the bandwidth in their sce- 
nario, which is an essential part of current OFDM standards. This as- 
sumption, causes the matrix Fp,cp to contradict the Restricted Iso- 
metric Property (RIP) defined in [?] and thus the use of Compressive 
Sensing (CS) algorithms as described in [?], unpractical. In the pres- 
ence of zero-padding, the matrix Fp,cp becomes ill-conditioned (a 
small variation in Hp leads to a great change in hcp). Due to zero- 
padding we do not have any pilots in the zero padded parts, further 
complicating the use of time-domain techniques. 



3. PROPOSED METHOD 

As previously stated, the OFDM Channel Impulse Response (CIR) is 
sparse. In this section we will propose a new scheme called Adaptive 
Thresholding for Sparse Signal Detection (ATSSD) that can exploit 
this inherent sparsity even with ill-conditioned matrices, Fp,cp- 

In this method, the spectrum of the channel is initially esti- 
mated using linear interpolation between pilot subcarriers. The 
zero-padded parts of the bandwidth which are possibly measured as 
nonzero due to the noise, are initially estimated as zero while at the 
middle subcarriers which are equipped with comb-type pilots, the 
mentioned linear interpolation using the noisy samples of the chan- 
nel at pilot locations is employed. This initial estimate is improved 
in a series of iterations that finds the sparsest time domain response. 

The initially estimated spectrum is passed through the IFFT 
block to obtain a crude version of the time domain impulse re- 
sponse. This initial estimate is fed to the ATSSD block which 
consists of several iterations. In each iteration this method tries 
to find the location of the taps via a thresholding scheme on the 
estimated channel from the previous iteration (or in the case of the 
first iteration, the crude version discussed above). The thresholding 
scheme is further described in sec. 13.11 After finding the location 
of the taps in each iteration (locations whose previously estimated 
amplitudes stay above the threshold), their respective amplitudes 
are again found using the Minimum Mean Square Error (MMSE) 
criterion discussed in sec. 13.21 In each iteration, due to thresholding, 
some of the fake taps which are noise samples whose amplitudes 
were above the threshold in the previous iteration, are discarded. 
Thus, the new iteration starts with a lower number of fake taps. 
Moreover, because of the MMSE estimator, the valid taps approach 
their actual values in each new iteration. In the last iteration, the ac- 
tual taps are detected and the MMSE estimator gives their respective 
values. The main steps of the proposed algorithm are: 

1. Using linear interpolation, set the crude time-domain channel 
as described above for the initial estimate. 

2. Discard the taps with amplitudes bellow the threshold. 

3. Estimate the value of the remaining taps using the MMSE 
criterion. 

4. Stop if the estimated channel is the same as the previous it- 
eration or when a maximum number of iterations is reached, 
else go to step 2. 

3.1. Channel Tap Detection via Thresholding 

As mentioned earlier, in each iteration, taps bellow a certain thresh- 
old are discarded. The threshold in the i*'' iteration is set as: 

Threshold = /3e°' (4) 

where i represents the iteration number; i.e., the threshold exponen- 
tially increases as the iterations proceed, a and /? are constants that 
depend on the number of taps and initial powers of noise and channel 
taps; since these parameters are hardly known, a rough approxima- 
tion is used. In the first iteration, the threshold is a small number and 
with each iteration it is gradually increased. Intuitively, this grad- 
ual increase of the threshold, results in a gradual reduction of fake 
taps (taps that are created due to noise). This is further explained in 
section|4]from a mathematical point of view. 



3.2. Estimation of Channel Tap Values via MMSE 

After finding time-domain channel tap positions, in eacli iteration we 
need to find an estimate of tlieir value. This is equivalent to solving 
a linear equation in the presence of noise. That is, we wish to obtain 
the value of the CIR vector h at tap positions (T), from the equation 



hT + W„ 



(5) 



where the vector Hp is the measured CFR vector at pilot positions, 
Fp,T is the sub-matrix of the DFT matrix (F) obtained by selecting 
the rows and columns that pertain to the pilot positions (p) and likely 
channel taps (T), respectively, and Wp is the noise vector at pilot 
positions. The MMSE estimator tries to find the estimate vector hr 
that minimizes i5{||hT — hT||2}, hence the name Minimum Mean 
Squared Error. This estimate vector is given by: 



iiT = Rh^F^y (Fp.TR-hrF^y + Rw) 



(6) 



where denotes the Hermitian operation and Rh^ and Rw are 
the auto-covariance matrices of hr and Wp respectively. When the 
noise vector W is a random complex white Gaussian process, Rw 
is equal to 2a'^I, where variance of both real and imaginary parts 
are assumed to be CTw- Also Rh^ can be approximated with Phi 
where Ph is the average power of the channel. Thus, the estimate 
vector can be written as: 



(7) 



where the channel power Ph can easily be estimated by dividing the 
power of the to-be-equalized data by the constellation power (power 
of the equalized data is the same as the average constellation power). 

4. MATHEMATICAL ANALYSIS 

In this section we present a mathematical intuition for the proposed 
method. For simplicity of the analysis, we assume that both channel 
taps and the additive noise in the time domain at the i*'' iteration (in- 
cluding the initial iteration formed by linear interpolation) are zero- 
mean normal complex random variables with equal real and imagi- 
nary variances tj^^jp and cr^^i, respectively. Therefore, their ampli- 
tudes which are of main concern in ATSSD have Rayleigh probabil- 
ity distribution function. In addition, we assume that the probability 
of a time sample being a channel tap is ptap', i.e., if Ncp represents 
the length of the cyclic prefix, we expect to have a channel with 
Ptap ■ Ncp taps. Now, if we set a threshold (rj) to distinguish the 
noise and valid samples, we probably miss some of the original taps 
while detecting a number of fake locations. The probability of False 
Alarm (FA); i.e., detection of fake taps,using the Rayleigh distribu- 
tion is given by: 



Pfa 



■ dx = e 



(8) 



Moreover, the average power of these fake taps (Etnp|„) is: 

'^tap\n = E{\wf I \w\ > 77} 
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Assuming 77 as the threshold for the i* iteration, we expect the over- 
all noise power (real and imaginary) of the next iteration to be equal 
to Etnp|„. Although this is only a rough statistical approximation, 
simulation results confirm that after MMSE estimation, the noise in 
the remaining taps (locations that stayed above the threshold in the 
previous iteration) is almost normal with zero mean and variance 
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For the choice of rj, in [?] it is shown that the threshold which maxi- 
mizes SNR of the signal after thresholding (considering both missed 
and fake taps) is found as: 
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Using ifTTTlinlfTOl) we get: 
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(12) 



which shows the increase in the noise variance; the noise samples 
which stay above the threshold should have higher amplitudes and 

therefore, higher variance. For reasonable SNR values of the chan- 

^2 

nel, — ^ is large enough to have the following approximations: 
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As can be seen in J13b . o-'^^i+i < ka'^ i\ however, we assume 
the upper bound for the noise and thus: 
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where /3 is a constant. Similar to l|4j, il4l shows the same exponen- 
tial increase in the threshold with respect to the iteration number. 

5. SIMULATION RESULTS 

To verify the efficacy of the proposed channel estimation, we per- 
formed computer simulations based on the DVB-H standard using 
MATLAB. The main simulation parameters and options are shown 
in Table [T] and |2] To show the robustness of ATSSD in estimating 
the CFR for time-invariant channels, we compared the BER after 
Viterbi decoding of the proposed method to that of the hypotheti- 
cal ideal channel in Fig. [TJ by ideal channel we mean that for data 
equalization, we used the exact channel instead of its estimate. As 
can be seen in Fig. [T] the BER performance of the proposed method 
almost coincides with that of the ideal channel. Thus, in this sense, 
the ATSSD estimation is perfect. Moreover, the BER performance 
of Linear Interpolation (a conventional channel estimation method) 



Table 1. Simulation Parameters 



Parameter 


Specifications 


DVB-H mode 


2K 


Number of carriers 


1705 


Number of pilot carriers 


146 


OFDM symbol duration 


224^s 


Gaurd Interval 


1/8 (256) 


Signal Constellation 


QAM- 16 


Convolutional Encoder 


niocT and 133oct 


Coding rate 


1/2 


[itermax , a , /3] 


[5 , 0.8 , 0.008] 


Table 2. Multi-Path Profile (Brazil Channel D) 



Delay (/is) 


0.0 


0.48 2.07 2.90 
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5.78 
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-0.1 
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Fig. 2. Performance of the proposed method at different Doppler 
frequencies for Brazil D channel. 



the same as if we used the exact channel even for time-varying chan- 
is also shown for comparison. The proposed method is also effec- nels. 
tive in time- varying channels. This fact is verified in Fig. |2] where hamash alakie 

the BER performance (after Viterbi decoding) of ATSSD is shown 
for different Doppler frequencies. As obvious in Fig. |2l ATSSD 
shows little performance degradation with increase of the Doppler 
frequency. 
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Fig. 1. Performance of the ideal, linear-interpolated and the pro- 
posed channel estimates under time-invariant Brazil D channel. 



6. CONCLUSION 

Due to the partial use of bandwidth in OFDM transmission, the re- 
duced frequency to time transformation is no longer stable; there- 
fore, most of the sparsity-based channel estimation methods diverge. 
We have proposed a thresholding method which solves this problem 
while benefiting from the sparsity criterion. Simulation results show 
that the performance of this method after data equalization is almost 



